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ABSTRACT 

We use a high-resolution A^-body simulation to investigate the influence of background galaxy 
properties, including redshift, size, shape and clustering, on the efficiency of forming giant arcs by 
gravitational lensing of rich galaxy clusters. Two large sets of ray-tracing simulations are carried out 
for 10 massive clusters at two redshifts, i.e. zi ~ 0.2 and 0.3. The virial mass (Mvir) of the simulated 
lens clusters at z ~ 0.2 ranges from 6.8 x lQ^^h~^MQ to 1.1 x IQ^^H'^Mq. The information of 
background galaxies brighter than 25 magnitude in the /-band is taken from Cosmological Evolution 
Survey (COSMOS) imaging data. Around 1.7 x 10^ strong lensing realizations with these images as 
background galaxies have been performed for each set. We find that the efficiency for forming giant 
arcs for z\ = 0.2 clusters is broadly consistent with observations. Our study on control source samples 
shows that the number of giant arcs is decreased by a factor of 1.05 and 1.61 when the COSMOS 
redshift distribution of galaxies is adopted, compared to the cases where all the galaxies were assumed 
to be in a single source plane at z\ = 1.0 and z\ = 1.5, respectively. We find that the efficiency 
of producing giant arcs by rich clusters is weakly dependent on the source size and clustering. Our 
principal finding is that a small proportion (~ 1/3) of galaxies with elongated shapes (e.g. ellipticity 
e = 1 — h/ a > 0.5) can boost the number of giant arcs substantially. Compared with recent studies 
where a uniform ellipticity distribution from to 0.5 is used for the sources, the adoption of directly 
observed shape distribution increases the number of giant arcs by a factor of ~ 2. Our results indicate 
that it is necessary to account for source information and survey parameters (such as point-spread- 
function, seeing) to make correct predictions of giant arcs and further to constrain the cosmological 
parameters. 

Subject headings: gravitational lensing - galaxies: clusters: general - dark matter - methods: data 
analysis 



1. INTRODUCTION 

Giant arcs are spectacular examples of strong gravitational lensing found in rich galaxy clusters. Background 
galaxies are stretched into long, thin arcs by the intense foreground gravitational field. Hundreds of giant arcs have 
been found in both optically-selected a nd X-ray selected clusters (e.g. iLuppino et al.l[T999tlZaritskv fc Gonzalez |[2003t 
IGladders et alllM)! ISa^d et al.l[2005h . 

Massive clusters are efficient producer of giant arcs and the number of giant arcs is a good indicator of the abun- 
dance of massive clusters. The halo mass function is very sensitive to the cosmological parameters, especially at the 
massive end. Moreover, the internal structures of massive clusters, such as substructures and ellipticity, also depend 

on the cos mologica l parameters. Th ey all affect the le nsing probability (optical depth, e.g. , ^eneghctti et al 2003aib| : 

Wambsaa nss et al.ir2 004: To rri et al.ll200l iDalal et"all[200liri et al.l[200alOguri et al.|[200 3. 2008; Hilb ert et al]l2007l : 
Puchwein et al !l !2005., 200ft) in various degree. Therefore, the observation o f giant a,rcs is a useful probe of the cosmo- 
logical model, in particular the matter power-spectrum normalization, erg (|Li et al.ll2006bf) . the matter density, rim,o, 
and to a less extent the cosmological constant, fl\_o. 

However, in order to use the observations of giant arcs to constrain the cosmological models, one has to thoroughly 
understand how the lensing probability is affected by the distribution of background source galaxies, as well as the 
intrinsic properties of lens population. It has been shown that the lensing probability increases significantly with 
the increase of the source redshift, thu s it is necessary to quantify the redshift distribution of source galaxies in 
lensing studi es ([Wa mbsgans s et al]|2004[) . The lensing efficiency of massive clusters does not depend on the source size 
si gnificantlv (iLi et a l. 2005), although not all the real source information has been used, in particular their shapes. 

iHoresh et al.l ()2d05l) first adopted realistic galaxy images with known photometric redshifts from the Hubble Deep 
Field (HDF) as background sources. They selected clusters from a cosmological A^-body simulation in a ~ 140/i~^Mpc 
box at zi ~ 0.2. The mass of the simulated lens clusters (0.54-1.1 x lO^^^h^^MQ) is similar to that of 10 X-ray-selected 
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clusters ()Sniith et al.ll2005f) . the (mean) mass range of which is around 6.3 x IO^^H^^Mq < M2Q < 2.0 x 10^^h~^Mc^ 
Note that the mass is calculated based on the X-ray luminosity by using the Lx-Af200 relation (jReiprich fc Bohringej 
120021) . They argued that the probability of producing giant arcs is 1 per cluster after observational effects are 
included, and emphasize that the number of giant arcs produced by the simulated clusters in their adopted cosmology 
(with ^ai,o — 0.3, r^A.o = 0-7 ^-nd the power-spectrum normalization of ag = 0.9) is consistent with the one observed 
in the massive clusters at redsliift 0.171 < Zc < 0.255. 

In this paper, we focus on the impact of background galaxy properties, rather than the properties of the lens clusters, 
on the efficiency of producing giant arcs. We use the COSMOS data as direct input for the properties of background 
galaxies. We investigate how the strong lensing probabilities are affected by the shape, size, r edshift distribution s 
and clustering of background galaxies. Instead of adopting all these source properties as a whole (i Horesh et al.l[2005l ). 
we vary each of the source properties in turn to gain a clear understanding about their individual impacts on strong 
lensing statistics. The COSMOS data has a substantial advantage over the HDF as it covers a much large r sky area 
(its sky coverage is 1.6 square degrees vs. 5.3 square arcmins for the HDF) as used in iHoresh et al.l (|2005D . and thus 
is less affected by the cosmic variance. While we confirm the effects of source redshift and size distributions, we 
find that the shape distribution of real galaxies has a dramatic effect on the lensing probability. With the COSMOS 
images, the lensing probability is boosted by a factor of 2 relative to that based on simple assumptions of the shape 
distribution (e.g. random ellipticity between and 0.5). Our results indicate that the shape of galaxies is an important 
factor in matching the theoretical predictions and the observations of giant arcs. We will also show that clustering of 
background galaxies has a negligible effect on the lensing efficiency of producing giant arcs. 

This paper is structured as follows. In §2, we discuss COSMOS and background source population. In §3, we discuss 
the simulation we use and our lensing methodology. In §4, we present our main results and compare with recent 
studies. We finish in §5 with a short summary and a brief discussion of the implications of our results. 

2. COSMOS AND BACKGROUND SOURCE POPULATION 

COSMOS is a deep survey covering a 2 square degrees of equatorial field containing over 2 million galaxies. It is an 
ideal sample of background galaxies for our lensing simulations due to the excellent resolution and large survey area (to 
reduce the cosmic variance). The image is centered on RA=10:00:28.6 and D EC=-|-02:12:21.0 from the Hubble Space 
Telescope Advanced Camera for Surveys (HST/ACS) (IK^ moer et al]l2007l ). Its calibration is relatively reliable due 
to the absence of atmospheric absorption. The /-band image data we use is the second public release of COSMOS 
observations. The camera has two 2048 x 4096 CCD chips with a pixel scale of 0.05". All the drizzled image data 
have been flux calibrated and astrometrically registered. The whole image in the /-band has been cut into 575 edge- 
overlapped tiles of 4960 x 4960 pixels of CCD, which are generated by rotating with an angle of ^ 10 degrees to the 
tile edge and embedded at the center of a rectangular box of 5600 x 5600 pixels. The left area in the box is empty of 
galaxies. 

The galaxy images in the COSMOS field are identified with the SExtractor program (jBertin fc Arnout^ Il996f l. 
The SExtractor is a well developed software program for extracting galaxy images from raw data and it produces a 
catalog of objects including properties such as their morphology, positions and magnitudes. The SExtractor detection 
threshold for all the ACS image data is set to be 1.5cr which is about 25mag/arcsec^ in surface brightness. Finally, 
a total of ~ 3 X 10^ galaxies are identified. The output of SExtractor contains the relative positions of sources in 
the tile and positions in the celestial equator coordinate system. The relative position is useful for identifying which 
galaxies are strongly lensed in the lensing process (se e ij3.2 | ). The ce lestial coordinates of galaxies are matched with 
the positions in COSMOS photometric catalog (Mobasher et al.ll2007f) of galaxies brighter than /-band AB magnitude 
of 25 to obtain the redshift information. Only the galaxies successfully matched, thus having correct photometric 
redshifts, are added in our lensing simulations. Using the COSMOS mask catalog, we exclude all star-like sources, 
galaxies close to star-like sources, and those possibly contaminated by spikes due to the overflow of CCD pixels. In 
particular, we use the same shape analysis program as in search for giant-arcs (see ^'6.2\i to identify very elongated 
sources with the major-axis to minor-axis ratio > 10 in all the tiles and exclude them to avoid "giant-arc like" artifact 
sources (e.g. imperfectly masked spikes due to CCD pixel overflow and asteroid trails) or real sources which could 
have been lensed by foreground objects in COSMOS field. 

All the cleaned COSMOS tiles are added behind the lens clusters as source patches. The sky coverage of strong 
lensing simulation is the same as the original COSMOS field. Due to the magnitude limit of the photometric catalog, 
the density of matched galaxies with photometric redshift is about 500 per tile of 5600 x 5600 pixels (~ 4.13 x 4.13 
square arcmins) which contains 4096 x 4096 effective CCD pixels. Thus, the surface number density of the extracted 
galaxies is about 43 galaxies arcmin"^ in the CCD imaging area. Tiles of 5600 x 5600 pixels with white margin are 
placed randomly behind the lensing clusters in the lensing simulation, without specifically requiring that the lensing 
caustic area falls within the actual image. This reduces the effective source number density to 23 galaxies arcmin"^ 
in our lensing analysis. 

The scatter plots of source size (in pixels) vs. redshift, the ellipticity (e) vs. redshift and ellipticity vs. source size 
with the median (as red filled circle) and the 68% value range (as red bars) for 20 randomly selected source tiles are 
shown in Figs. [T][3l The ellipticity e is quantified by 1 — b/a, where b and a are semi-minor and semi-major axes of a 
galaxy image. The axes are obtained by fitting the image with an ellipse using the same method as that for quantifying 
the shape of giant arcs (see ij3.2p . Fig.[l]shows that the sources become smaller as the redshift increases from to ~ 2, 
and then remain flat out to redshift z = 3.0. However, the ellipticity of galaxies is almost independent of redshift (see 
Fig. [2]). The source size and ellipticity are somewhat correlated (especially for the sources of the effective diameter 
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smaller than ~ 1.2", as can be seen in Fig. larger sources appear to be somewhat more elliptical. 

The photometric redshift distribution of COSMOS background galaxies is plotted (as the black solid line) in Fig. [4] 
for all the galaxies in 575 tiles of COSMOS field. As can be seen the redshift distribution peaks around z = 0.8, but 
extends out to z = 3.0 with a median redshift of 0.78. The histogram is well fitted by the empirical formula from 
ISmail et al.l (Il995l) : 



/3 



P^^^^ " 3T^/Q/m exp - — , (1) 

where zq ^ 0.5 and /3 ~ 1.2; the best- fit line is indicated as the black dashed line in Fig. H) 

With the photometric redshifts, we can construct an approximately three-dimensional distribution of background 
galaxies in order to study the effect of source redshift distribution on strong lensing efficiency. A series of source planes 
perpendicular to the line of sight is set between redshift z\ ~ 0.2 0.3) to Zg — 3.0. The redshift separation between 
two adjacent source planes dz is set to 6z/z = 0.1 for 1.5 < z < 3.0 and to 6z/z = 0.05 for 0.2 < z < 1.5. The galaxies 
are projected onto the planes according to which redshift bin they fall in. In each source plane, we mark all the galaxy 
pixels, which are used to identify the lensed sources in subsequent ray-tracing. There are 48 source planes in total for 
the ray-tracing program (see §3.2p . 

We will also study the impact of ellipticity and size distributions on lensing cross-sections. The ellipticity distribution 
of the COSMOS galaxies are illustrated (as the black s olid histogram) i n Fig . [5l The ellipticity I ~ b/a distribution 
has a median of 0.4, which is similar to that presented in [Ferguson et al] ()2004[ ). Compared to the studies of giant arcs 



using a uniform distribution of the ellipticity with < 1 — b/a < 0.5 for background sources, COSMOS clearly has an 
extended tail of highly flattened galaxies, which will have a significant impact on the efficiency of forming giant arcs 
(see i |4X2l) . 

The size of each galaxy is quantified by an effective diameter Doff defined as the diameter of a circular-shaped galaxy 
with the same number of the image pixels. The size distribution of galaxies in these tiles is shown (as the black 
solid line) in Fig. [6] The size distribution of galaxies has a peak at D^s ~ 0.75", but ranges from effective diameter 
Deff < 0.2" t o Dpff > 2". T he median diameter of background galaxies is about 0.6", smaller than that assumed to 
be Doff = 1" (|Li et al.l 120051 ) since a relatively deeper galaxy survey is adopted in this work. 

3. NUMERICAL SIMULATION AND LENSING METHODOLOGY 

3.1. Simulated Clusters 

The l ens clusters a r e selec t ed from a A^-body simul ation generated with a upgraded version of vectorized-parallel P'^M 
code of iJing k Sutol (|2002[) (|Jing. Suto k Mol 120071 ). The underlying cosmological model is a ACDM (rj^^o = 0.268, 
^^A,o = 0.732, n = 1) model. The simulation uses 1024^^ dark matter particles in a box with a side length of 600/i~^Mpc 
(comoving). The mass resolution is 1.5 x Mq, thus m assive clusters are reason ably resolved (with more than 

4 X 10^ particles in our samples). In this simulation (refer to lJing. Suto fc Mol (2007) for more details), the Hubble 
constant h in unit of 100 km s^^ Mpc^^ is take n to be 0.71, while the ampl itude of the linear density power spectrum, 
as, is 0.85, which is obtained from CMBFAST (jSeliak fc Zaldarriagalll996l ) directly. Particle pairwise interactions are 
softened on scales smaller than 30 /i~^kpc. The clusters are identified by the Friends-Of-Friends (FOF) method with 
a linking length of 0.2 times the mean particle separation, and the cluster mass Myjr is defined as the mass enclosed 
within the virial radius according to the spherical collapse model ()Kitavama fc Sutol IT996I: IBrvan fc Nof^m^ HOOSl : 
I.Ting fc Sutol [20(il ). 

We primarily use a sample of 10 massive clusters at redshift z ~ 0.2 to perf orm ray-tracing si mulations. This is 
to facilitate comparisons with observations of the X-ray selected cluster sample ([Smith et al.l 120051 ). The mass of the 
massive clusters at redshift z ^ 0.2 ranges from 6.8 x 10^*h~^MQ to 1.1 x 10^^h~^ Mq, roughly consistent with the 
observed range of mass (6.3 x lO" < M200 < 2.0 x IO^^Ii-^Mq of clusters at redshift 0.171 < Zc < 0.255).^ We shall 
also use 10 most massive clusters at z ~ 0.3 for lensing efficiency comparisons. 

3.2. Ray-tracing and Mapping 

We use the thin lens plane approximation for t he lensing calculation; the line of sight contributions may be important 
(jWambsganss et all 120041 : iPuchwein et al.ir2009f ). The surface density is calculated on a mesh of 102 4 x 1024 (coarse) 
grids centered at the cluster with a size of 2 times the virial radii. Our sm oothing method, similar to lLi et all (12006a), 
uses the smoothed particle hydrodynamics kernel (jMonaghan et al.]|1992f ) to distribute the mass of particles enclosed 
inside the virial radius of the lens onto 3-D grids and then the surface density is obtained by integrating along the line 
of sight. The dimensionless surface density (k) map is obtained by dividing the surface density by the critical value, 

D 

(2) 



47rG DiDis ' 

where Dg, Di, and Dig are the angular diameter distances between the source and the observer, the lens and the 
observer, and the source and the lens. The lensing potential on the grids of the coarse mesh, 0, is calculated by the 

^ Notice that M200 is defined as the mass enclosed within in r200 1 the radius within which the average density is equal to 200 times the 
critical cosmological density. The difference between M200 and A/vir are within 10% for 8 of 10 clusters. It is within 30% for the other two 
clusters. 
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Redshift 

Fig. 1. — : The scatter plot of size against redshift for galaxies randomly selected from 20 COSMOS image tiles with 
a bin width of Az = 0.3. Each tile has 4096 x 4096 effective pixels (with a pixel scale of 0.05 arcsec by 0.05 arcsec) 
covering a field of view of 3.4 arcmins by 3.4 arcmins. The median (red filled circle) and the 68% source size range 
(red bar) are indicated in 10 b ins of redshift. As can be seen, the source size distribution is almost independent of the 
redshift. For comparison with [Ferguson et al.l (|2004[ ). the mean value (blue filled circle) of effective diameter in each 
redshift bin is also plotted. As our survey depth is much shallower, they actually agree well with each other. 

Fast Fourier Transform (FFT) method from the dimensionless surface density map since they are related by V^(/) = 2k. 
The shear (7) is then obtained by the second-order derivatives of the potential, and the (signed) magnification is given 
by M-i = (1-^)2-72. 

To better determine the geometry of an image and to perform an efficient lensing simulation, we refine the cell size of 
the central high magnification area with > 2.5 in the lens plane up to 0.05", identical to the intrinsic observational 
resolution in the source planes. High resolution is necessary for resolving the lensed image well, especially in the 
direction of minor axis, which is much more sensitive to the resolution. The lensing potential on the refined mesh is 
obtained by cubic spline interpolations of the neighboring 14 x 14 coarse grids in the 1024 x 1024 mesh. Once we 
obtain the lensing potential in the whole fine mesh area of interests, the deflection angle on each grid point can be 
calculated from the derivative of the lensing potential, a — V0. The lens mapping from the lens plane to each source 
plane at different redshifts can be constructed using the lens equation (i.e. ray-tracing), y = x — a, where y, x are 
dimensionless source position in source plane and image position in lens plane. The ray-tracing procedure is repeated 
for each of the ~ 50 source planes, as we discussed in f}2] 

We identify the image pixels in the lens plane from the mapping between the lens plane and the source planes. 
For each projection of a lens, we go through all the refined grid points one by one and check whether there is a 
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Redshift 

Fig. 2. — : The scatter plot of the ellipticity, e = 1 — 6/a, against redshift for galaxies randomly selected from 20 tiles 
with a bin width of Az = 0.3. The median (red filled circle) and the 68% ellipticity range (red bar) are indicated in 
10 bins of redshift. As can be seen, the source ellipticity distribution is randomly distributed at all the redshifts. The 
mean value (blue filled circle) of ellipticity in each redshift bin is also plotted for comparison. 

corresponding point in the source planes falling into a galaxy. If so, we identify this grid point in the image plane 
as a lensed image pixel of that source galaxy; if not, it is a trivial grid. All the highly magnified lensed image pixels 
(e.g., > 2.5) can be located. Finally, we can obtain the whole image(s) of a background galaxy by collecting all the 
image pixels of the same source using Friends-of-Friends (FOF) method. 

For source slices at different redshifts, the high-magnification areas will be mapped into areas of different sizes in 
the sky. We choose a rectangular conjunctive angular area inclusive of every caustic area and put this angular box 
randomly in the source tile to select different background galaxies for the ray-tracing simulation. The box size is 
usually on the level of a few square arcmins (depending on projection and also mass profile). Since each added source 
tile is about 4.7' x 4.7', around 10 random lensing realizations are needed to cover a tile effectively. To avoid boundary 
effects, the source tile is periodically expanded in both dimensions which are then randomly sampled. 

The same process is carried out for all the three orthogonal projections of the massive clusters at the expected 
redshifts. 10 ray-tracing simulations are carried out with all source galaxy samples in the COSMOS field for each 
of 3 orthogonal projections. The lensing images found are stored and subject to more detailed analysis to find their 
length-to- width ( L/W) ratios, sizes, the r edshi fts. Of particu lar importance is the L/W ratio. We determine this 
quantity following iBartelmann et all (|1998 ^ and 'L i et al.l (|2005f ). Briefly, we first identify the center of the image pixels 
for a given source (xc), and then locate the point xi that is the furthest from the centre, and finally the point X2 that 
is the furthest from xi. We then fit a circular arc that passes through these three points (xi, X2, Xc). The half length 
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Effective Diameter D^ff/arcsec 

Fig. 3. — : The scatter plot of the eUipticity, t — 1 — 6/a, against size of galaxies in 20 randomly selected tiles with 
a bin width of ADoff = 0.3. The median (red filled circle) and the 68% eUipticity range (red bars) are indicated in 
10 bins of effective diameter Deff. The galaxies are more elliptical as the size increases for sources with the effective 
diameter smaller than ~ 1.2", while the median eUipticity keeps around 0.5 for the left size range. The mean value 
(blue filled circle) of eUipticity in each effective diameter bin is also plotted for comparison. 

of the arc is taken to be the semi-major axis length of the ellipse, a, and the length of the semi-minor axis is taken to 
be 6 = S'imago/(7i'a), where 5'imago is the area covered by the image. The length-to-width ratio is then simply ajh. As 
mentioned before, the same procedure is used to determine the eUipticity, 1 — 6/a, for unlensed galaxies. 

3.3. Simulating test images 

We carried out a test by putting a source around the tangential caustic curve (as the green dotted line) of a massive 
simulated cluster. The simulated images (shown in red) are illustrated in Fig. [7] together with the source (shown in 
black) positions. As expected, it produces one, two and three images when the mocked source is outside, on and inside 
the caustic curve. Fig. [8] shows a simulated image of a massive lens cluster of a realization based on a real source tile. 
Clearly, the sources (shown in black) have been systematically stretched around the cluster center. In particular, there 
is one giant arc at redshift z = 2.52 (labeled as "Al") with L/W ^ 22 among the lensed images produced from source 
"a" , which has an additional image labeled as "A2" . Notice that the positions of the images and background galaxies 
in different planes are transformed into the same coordinate system centered at the cluster centre. 

4. RESULTS 

Many recent studies have shown that the efficiency of forming giant arcs depends on the background source popu- 
lation, in particular, their sizes and redshifts fe.g.. .Wambsganss et al..,2004; ,Li et al...2005; .Horesh et al.,,2005i) . Since 
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Fig. 4. — : Comparison plot of redshift probability distributions between the subpopulation of original galaxies which 
have been strongly lensed to be arcs of L/W> 10 {zi = 0.2) and the parent population of COSMOS galaxies. The 
number fractions per redshift bin of width Az — 0.2 are plotted. The lensed population is clearly shifted toward 
high-redshift end. The fitting curve of the COSMOS redshift by Eq. [T] is indicated as the black dashed line. The 
redshift distribution of the fiducial control sample consisting of 100 randomly selected tiles (used in Table [1]), which 
roughly follows the mea n di stribution, is plott ed in magenta. For comparison, the redshift distributions used in 
iBartelmann et all (|1998f l and lFedeh etall (|2008[ ) are also shown. 

different giant arc surveys have different selection effects, the goal of this work is not only to compare our strong 
lensing efficiency with observations, but also to focus on using realistic, empirical information of background sources, 
and then to study how the lensing efficiency changes as we vary the assumptions. 

For this purpose, we construct 7 comparison source samples in the control group which will be added in the ray- 
tracing simulations. All the samples are based on 100 randomly selected COSMOS tiles. The mean number density 
of this control sample is around 500 galaxies per tile, which is the same as the mean density from the whole field. The 
sample with COSMOS background sources is our fiducial case (Case 1). This sample also provides the basis on which 
we generate mocked samples with partially or fully artificial source properties for Cases from 2 to 7. The artificial 
sources are mocked tile by tile. The lensing simulations are performed for 3 projections of 5 massive clusters. For 
effectively sampling the source area, 10 realizations are carried out for each projection in each source tile. There are 
15000 simulations in total for each case. The resulting arc statistics for this control group (e.g. six mocked samples 
together with the COSMOS sample) are listed in Table [TJ The lensing efficiency changes as different source properties 
are switched. We find source properties affects the lensing efficiency in various degrees by comparing the numbers of 
giant arcs between different cases, especially the source ellipticity and redshift. Furthermore, we also find that the 
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Fig. 5. — : Comparison plot of shape probability distributions between the subpopulation of original galaxies which 
have been strongly lensed to be arcs of L/W> 10 {zi = 0.2) and the parent population of COSMOS galaxies. The 
number fractions per ellipticity bin of width Ae = 0.1 are plotted. The lensed population clearly shows a bias towards 
high ellipticity compared with the original distribution. The shape distribution of the fiducial control sample consisting 
of 100 randomly selected tiles (used in Table [H) is pl o tted i n magenta and agrees with COSMOS distribution very 
well. For comparison, the distribution used in lLi et all (|2005D is also shown. 

boosting factor due to the source ellipticity on lensing efficiency depends on the source redshift. The details are listed 
below. 



4.1. Lensing efficiency of COSMOS sources 

Using the empirical COSMOS background source information, and the method outlined in the last section, we 
implement a mock lensing survey for each projection of the massive clusters at reds hift zi = 0.2 and z\ — 0.3 for all 
of 575 tiles from COSMOS (Table [2). We adopt the definition of the L/W ratio of iLi et all (2005). In total, 12816 
giant arcs of L/W > 10 are found out of 172500 cases at z\ = 0.2, while 9704 giant arcs for clusters at z\ — 0.3. 
Therefore, the lensing efficiency is about 0.0743 giant arcs per cluster at z\ = 0.2 and 0.0563 giant arcs per cluster 
at zi = 0.3. We illustrate the lens redshift (zi) dependence and the fluctuation of number counting in Fig. [9] for 20 
randomly selected tiles. All the giant arcs {L/W > 10) produced by the 10 massive lens clusters from the galaxies in 
each tile are taken into account. The lensing probability appears to decrease from zi = 0.2 to zi = 0.5 in our case, 
although the fluctuation is relatively large mainly due to the difference of source properties between each tile and the 
intrinsic surface density fluctuations (e.g., between actual imaging region and white margin) within the tiles. 
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Effective Diameter D^jf/arcsec 

Fig. 6. — : Comparison plot of size probability distributions between the subpopulation of original galaxies which have 
been strongly lensed to be arcs of L/W> 10 (zi — 0.2) and parent population of COSMOS galaxies. The number 
fractions per effective diameter bin of width AD^s = 0.25 arcsec are plotted. There is little difference between the 
two distributions, which indicates they share the almost same population. The size distribution of the fiducial control 
sample consisting of 100 randomly selected tiles (used in Table [J) is plotted in magenta. There is no big d ifference 
from t he mean in most _Dcff range, except the places around the peak. For comparison, the distribution used in lLi et aTl 
(|2005f l is also indicated. 

4.1.1. Redshift influence 

In recent works (jBartelmann et aLlll998HLi et al.ll2005[ ). the background sources are taken to be at the same redshift, 
e.g., z = 1, to estimate the theoretical overall efficiency of forming giant arcs in our universe. As shown in Fig. 2] the 
background sources in fact have a rather broad distribution, which we take into account by binning them into ^ 50 
redshift slices. If we put all the sources into a single redshift plane at Zg = 1, we find that in this case the number 
of giant arcs {L/W > 10) increases by a factor of ~ 1.05 compared with Case 1 in Table [l] since the equivalent 
source redshift plane, in which we put all the source galaxies to generate similar number of giant arcs, for the 100 
randomly selected COSMOS tiles (i.e., "fiducial control sample" tiles which are considered as the "seeds" to generate 
mock samples for cases in Table [T]) is about 0.9 which is slightly lower than Zs = 1.0. Notice that the shape, size 
distributions and surface density follow the original distributions of COSMOS. Therefor e, Zg pla ne is coinci dently a 
good approximation to the broad redshift distribution. However, according to Wambsgan ss et all |2004) and lLi et al.l 
|2005), the lensing probability becomes larger when the sources are put at a higher redshift, which is also seen in the 
redshift distribution of lensed sources in Fig. [4] where the lensed galaxies, stretched to become giant arcs oi L/W > 10, 
are clearly biased toward high redshift relative to the underlying background population. For comparison, the analytical 
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Fig. 7. — : Simulated images for a mocked circular galaxy located outside, on and inside the tangential caustic of 
a merging galaxy cluster, which is selected to illustrate for easy identification. The black dots indicate the sources, 
while the red ones are the lensed images. The green and blue dotted lines indicate the caustics and critical curves 
respectively. As can be seen, the lens produces single, double and triple strongly lensed images. 



TABLE 1: Results of 7 different sets of simulations, "z", "s"and "e" indicate the redshift , size and ellipticity of 
the background sources, while "l"and "0" represent following COSMOS and mocked. For example, "eO" implies the 
ellipticity (1-b/a) is randomly chosen in the range between to 0.5. We run 10 lensing simulations in each projection 
of the 5 (middle) massive clusters at zi = 0.2 for each source tile. In each case, 100 tiles are selected to construct the 
control sample with required source properties. The redshift, ellipticity and size distributions of the selected COSMOS 
tiles are plotted in Fig. [4l [5] and[6l labeled as "Control". There are 15000 ray-tracing realizations performed for each 
control sample. (Similar results are obtained for clusters at zi = 0.3). 



Case 


Name 


z 


s 


e 


N 


1 


zlslel 


COSMOS 


COSMOS 


COSMOS 


1204 


2 


zOslel 


z=1.0 


COSMOS 


COSMOS 


1264 


3 


zOsOeO 


z=1.0 


Dcff = 1-0" 


RANDOM 


321 


4 


zlsOel 


COSMOS 


D,a = 1.0" 


COSMOS 


852 


5 


zlsOeO 


COSMOS 


D,ft = 1.0" 


RANDOM 


369 


6 


zlsleO 


COSMOS 


COSMOS 


RANDOM 


493 


7 


zOslel 


z=1.5 


COSMOS 


COSMOS 


1935 



TABLE 2: Ray-tracing simulations are carried out for the massive clusters selected from two snapshots at redshift 
~ 0.2 and 0.3. The cosmological model is a ACDM model with fl^^ = 0.268, I^a.q = 0.732, n = 1, and as = 0.85. 
We implement 10 lensing simulations (marked as Nn) for each projection and each source tile. The whole simulation 
covers the full sample of 575 source tiles (marked as Nt). 



Set 




Mvir{lO^^/»-^M0) 


Nc 




Nr 


Nt 


1 


0.2 


6.8-11 


10 


3 


10 


575 


2 


0.3 


6.0-8.1 


10 


3 


10 


575 



redshift distribution in lFedeli et al.l (I2008D is a lso shown in Fig. HI If we adopt a source redshift distribution following a 
deeper galaxy survey like in lFedeli et al.l (|2008[ ) (keeping other properties unchanged) , the giant-arc number is increased 
by a factor ^ 2, compared with the number of giant arcs when the COSMOS redshift distribution is used (Case 1 in 
Table [J). 

However, the equivalent source redshift plane changes with source shape distributions. For example, considering the 
comparison cases. Cases 3 and 5 in Table [l] wc find a higher equivalent source redshift (e.g., Zs ~ 1.6) for galaxies of 
simulated shapes as in Case 3 than those of the intrinsic shape distribution in Case 1, which share the same intrinsic 
redshift distribution. This can also be explained by Fig. [TOl Galaxies of simulated ellipticities at higher redshift have 
a relatively higher weight than those of intrinsic ellipticities in producing giant arcs. It causes the equivalent source 
redshift for Case 5 to be larger than the mocked redshift in Case 3, e.g. Zs = 1.0, thus a higher lensing efficiency. 
Also, the redshift of equivalent source plane may shift to a higher redshift for a deeper survey. In this case, the lensing 
probability may increase substantially and a single source plane approximation (at Zg = 1.0) can become violated. For 
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Fig. 8. — : Simulated images by a massive lens cluster {z\ = 0.3) for galaxies in a COSMOS tile. The sources and their 
lensed images are shown in black and red respectively. The other unlensed galaxies in the tile are not plotted together 
for a clear view . The blue and green curves are critical lines and their corresponding caustics for sources at redshift 
Zs = 3.0, while the cyan cross indicates the center of lens. A strongly stretched arc "Al" is clearly seen corresponding 
to the source "a" ; an additional image "A2" is highlighted (in magenta) for the same source. 

example, we find an in crease of a factor ^ 2 in giant-arc number for a mocked source sample following the analytical 
redshift distribution inl Fedeli et al.l (|2008f ). which is shown in Fig. [H It is necessary to take account of the full redshift 
distribution in predicting the strong lensing efficiency. 



4.1.2. Source shape influence 

Previous studies (|Li et al.l 120051 : IPuchwein et al.l 120091 ) often assumed the sources follow a uniform distribution in 
the ellipticity (e = 1 — 6/a) between and 0.5. However, as can be seen in Fig. [5l while it is true that the ellipticity 
1 — b/a of most sources is less than 0.5, there is still a high fraction of sources, about 1/3, with ellipticity larger than 
0.5. The sources with large ellipticities, such as edge-on spiral galaxies, could be stretched much more easily to form 
giant arcs of large L/W values when the tangential direction of the lens is approximately along the direction of the 
major axis of the background galaxy. Therefore, the small portion of high ellipticity sources may change the giant-arc 
statistics substantially. 

For all the COSMOS sources in the fiducial control tiles, we change the shape of sources with random ellipticities (be- 
tween and 0.5) and orientations (between and tt) while keeping the source positions, sizes and redshifts unchanged. 
We perform a lensing simulation for this background galaxy sample. This case (Case 6) is labeled as "zlsleO". We 
find that the number of giant arcs is only ^ 1/2 of the value from the COSMOS sources. The shape influence is also 
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Tiles 

Fig. 9. — : The giant-arc statistics for three test runs at different lens redshifts, e.g., zi = 0.2, 0.3 and 0.5 are illustrated. 
The 20 test source tiles are randomly selected. For each source tile, all the giant arcs produced in 300 realizations of 
10 massive clusters are counted (see Table [5]). As can be seen, the fluctuation is large for a single source tile in the 
mock survey due to the limited realizations. The lensing efficiency decreases when the lens redshift ranges from 0.2 to 
0.5, which is consistent with the result in i )4.5l 

directly reflected in the different ellipticity distributions of the original COSMOS galaxies and the lensed population. 
As clearly seen in Fig. [5l the lensed population of original galaxies with giant arcs of L/W > 10 shows a clear shift to 
the higher ellipticity tail. 

It is important to explore whether the effect of source shape changes with the source redshift. Although the overall 
cross-section becomes larger for a higher source redshift, the change of the effective lensing area, in which elongated 
sources can be lensed to giant arcs while the simulated ones of < e < 0.5 cannot be, is not clear. We construct 4 
comparison pairs of control source samples generated from the fiducial control tiles, i.e. "zOsOel" and "zOsOeO", at 
redshift of Zg = 0.8, 1.0, 1.5, 3.5. As clearly seen in Fig. [TOl while the number of giant arcs increases with the source 
redshift, the shape impact on lensing efficiency decreases from ~ 4 to ~ 2, when source redshift ranges from Zs=0.8 
to 3.5. The boost effect of source shape is expected to be larger at Zg lower than 1.5 since the lensing cross-section is 
smaller at lower redshift and thus the enhancement due to the elliptical sources is relatively more important. 

4.1.3. Source size influence 

Instead of assuming a same effective diameter for all sources (jLi et al.ll2005l : IPuchwein et al.ll2009( ). we have adopted 
the COSMOS galaxy sizes in our fiducial case (Case 1). To see the size effect of background galaxies, we first check 
the size distribution for the lensed galaxies which are stretched to be giant arcs of L/W > 10. The comparison results 
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Fig. 10. — : The shape impact on cross-section dependence of Zs is illustrated. It is clearly seen that the lensing 
incidence is boosted after the real galaxy shape is adopted (top panel) . The Poisson fluctuations are indicated as error 
bars in the panel. As a comparison, we put the sources of random ellipticity (shown in red) or real ellipticity (shown 
in black) in source planes at redshift z = 0.8, z = 1.0, z = 1.5 and z = 3.5, keeping source size of Deff = 1". While 
confirming the shape impact on lensing cross-section is about a factor of 2 for Zs > 1.5, with a larger factor at z < 1 
(bottom panel). 

are illustrated in Fig. [6l it shows that there is little difference between these two distributions, and thus the source 
size is expected to have small effect on producing giant arcs. 

This is directly confirmed by two comparison samples, where we fix the effective diameters assigned to be Dcs ~ 1.0" 
and DcS = 0.5". We find a similar number of giant arcs (188 and 166 out of 3000 realizations), which is consistent 
with Fig. [S] Source size thus has a moderate effect on producing giant arcs when massive clusters are taken as lenses 
(also see Fig. [3] and Fig. [TT|) . The reason is that the typical source size is much smaller than the strong lensing caustic 
size, and thus a moderate change in the source size has a small impact on the number of giant arcs. 



4.2. Combination of Morphology and Redshift Distribution 

To see the overall impact of the simple assumptions for the shape, size, and redshift of background galaxies, we 
combine these effects to do lensing simulations (for random position effect, see > j4.3p . The redshift of all these sources 
is assigned to be 1.0. The shape with a random orientation (between and tt) are simulated as above, while the 
effective diameter is assigned to be 1". We find that the lensing efficiency for the Case 3 ("zOsOeO") is a factor of 3.75 
lower than that for Case 1 ("zlslel"). It is roughly consistent with the results from the combination of the individual 
effects, i.e. redshift effect (a factor of 1.15 which is got from comparing Case 3 and 5), shape effect (a factor of 2.31 
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Fig. 11. — : The color- weighted distribution plot of the ellipticity against size for sources lensed to be images with 
L/W > 10. The lensed sample is obtained from tracing back all the giant arcs produced by 10 massive clusters at 
redshift z ~ 0.2. The pattern of the distribution is similar to that shown in Fig. [31 although the elongated sources 
have a higher weight, which is consistent with that implied in Fig.O 

from Case 4 and 5) and size effect (a factor of 1.34 from Case 5 and 6), which is 1.15 x 2.31 x 1.34 = 3.56. It thus 
indicates that the galaxy shape, redshift and size are affecting strong lensing efficiency independently as well as these 
three properties are independent of each other intrinsically as shown in Fig. [1] [2] and [3] 

4.3. Effect of Source Clustering 

Theoretical calculations often assume that the background sources are randomly distributed in space. However, 
background galaxies are clearly clustered. It is therefore interesting to check the effect of source clustering on the 
incidence of giant arcs. 

The influence of the clustering effect is investigated by carrying out a ray-tracing simulation on a control sample of 
galaxies where the positions of original COSMOS galaxies in a fiducial tile are randomly shuffled while keeping the 
source surface density, size and ellipticity unchanged. We find that the efficiency of producing giant arcs (166 giant 
arcs out of 3000 realizations for 10 clusters at z\ = 0.3) is similar to what obtained in the COSMOS case (163 giant 
arcs) which implies the clustering effect of background sources is negligible on strong lensing. The likely reason is that 
the intrinsic source density is too low (43 per square arcmin) and close pairs of sources are too few to form connected 
giant arcs or to be strongly lensed simultaneously. Our results suggest that sources can be assumed to be random in 
the prediction of giant arcs, at least to the depth of COSMOS (in deeper images with a larger source density, clustering 
may become more important). 
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4.4. Effects of Seeing in Ground-based Observations 

For our fiducial simulation (Case 1), the observational data are taken by HST (as in the Hubble Deep Field used in 
Horesh et al. 2005). However, many giant arcs were discovered from the ground-based telescopes where the effect of 
seeing may be important. A proper evaluation would require a detailed simulation of seeing effects. Here we discuss 
briefly the approximate seeing effect on the predicted number of giant arcs. 

In COSMOS, there is a significant portion (~ 80%) of small galaxies with the effective diameter DeS < 1" (see Fig.|6]). 
The width of lensed giant arcs produced from these small sources probably remains less than l'^ In such cases ^ seeing 
(such as ~ 1") will have a much larger effect on the width of giant arcs than on the length (jMene ghetti et al. 2007), 
and in general the L/W ratio will be substantially reduced. Assume the seeing is ~ 0.5", the width of arcs intrinsically 
smaller than 0.5" will be blurred to ~ 0.5" which we adopt as the "observed" width. We find that in this case the 
fraction of giant arcs of L/W > 10 will be reduced to 80% of the COSMOS value. If seeing is as bad as 1", the lensing 
probability could be reduced to about 40% of the original value. Although this is only an approximate estimation, it is 
clear that the impact of seeing should be carefully evaluated in a detailed comparison between theoretical predictions 
and ground-based giant-arc surveys. 

4.5. Comparison with Observations and Horesh et al. (2005) 

The giant arcs in this work are selected b y quanti f ying t he ratio of length to width {L/W > 10). The definition 
of L/W in this paper is the same as that in Li et al. (2005). As a mock strong lensing survey, we have done 172500 
realizations of ray-tracing simulations for all the 10 lens clusters at z\ = 0.2 or zi = 0.3 with 575 galaxy tiles, and find 
that the lensing efficiency is about ^ 0.0743 arcs for clusters at z\ — 0.2 and ~ 0.0563 arcs for clusters at z\ — 0.3 per 
realization. 

On the face valu e, they are much lower than the expected mean probability of ^ 1 giant arc per realization 
(jHoresh et al.l[2005h . However, they adopted galaxy images in the Hubble Deep Field whose source density is about 
12 times higher (1378/5.3=260 arcmin"^ in HDF, vs. 23 arcmin^^ in our ca se). The lens i ng ef ficiency with the 
surface density normalized is 0.840 arcs and 0.636 per realization, respectively. iHoresh et al.l ()2005( ) used a different 
definition of L/W ratio which is smaller than ours by a factor of A/tt (e.g., their L/W > 10 is equivalent to our 
L/W > 2.57r ~ 7.5), and thus our lensing efficiency (of L/W > 10) will be boosted by a factor of (4/7r)^ if we adopt 
their length-to-width ratio definition. Taken these two factors into account, our clusters produce 0.840 x (4/7r)^ = 1.36 
and 0.636 x (4/7r)^ = 1.03 giant arcs per realization for z\ = 0.2 and 0.3, respectively. In this exercise, we have assumed 
that the giant arc formation cross-section scales as the magnification probability distribution {p{> \^\) cx |^| » 1). 
Therefore, the capability of our simulated lens clusters is in fact consistent with the observation, which is about 1.2 
arcs per cluster found by Horesh et al. (2005) using the sample in Smith et al. (2005). 

The slight excess in the number counting could be induced by small number statistics and the relatively narrow 
mass range of the simulated lenses by contrast to the observational clusters. It also could be due to the different lens 
redshift distributions of the simulated (zi ^ 0.2) and observational samples (0.171 < Zc < 0.255). 

5. SUMMARY AND DISCUSSION 

In this paper, aiming to study the impacts of background sources on strong lensing statistics, we use the /-band 
galaxy image data of HST/ ACS in the COSMOS to quantify the distributions of background source size, shape and 
redshift. Each galaxy image is extracted by SExtractor limiting a surface brightness down to ^ 25mag/arcsec^. The 
redshift of each galaxy image is obtained by matching its celestial position with the COSMOS photometric catalog. The 
selected sources are then lensed by 10 massive clusters of the mass range 6.8 x lO^'^h~^MQ < Mvir < 1-1 x lO^^h~^M0 
at zi ~ 0.2 and 6.0 x 10^'*/i~^Mq < Mvir < 8.1 x 10^'*h~^MQ at zi ~ 0.3 as the lensing clusters chosen from a 
cosmological simulation in the ACDM cosmology {flm,o = 0.268, r^A^o — 0.732). 575 source tiles within the COSMOS 
field of around 2 deg^ are fully used for the statistic study. 172500 ray-tracing lensing simulations are carried out for 
10 lens clusters to reduce the statistical fluctuation at the expected lens redshifts (i.e. zi = 0.2 and 0. 3 ). Th e incidence 
of giant-arc production in our simulation is roughly consistent with that observed in iSmith et al.l ()2005D , after the 
density difference is taken into account. The impacts of source size, shape and redshift on strong lensing statistics are 
investigated in detail. 

We find that the source size (less than a factor of 1.5) and clustering only have small effects on the production of 
giant arcs. In contrast, the dependencies on the source redshift and ellipticity are much more significant. The first was 
highlighted by Wambsganss et all (|2004f ). while the second is the main new finding of the current work. We find that 
adopting the empirical ellipticity distribution of COSMOS increases the lensing probability by a factor of 2 (see Fig.fTOl 
and also the number ratio between Case 1 and 6). The boosting effect of the ellipticity of the background galaxies has 
not been emphasized in the previous works. It should be included in theoretical modeling of giant arcs in future. This 
may also be of particular relevance for setting constraints on the power-spectrum normalization parameter as- 

There are a number of limitations in the present work. When scaling our results with the source number density, 
we implicitly assume that other source properties (like ellipticity and redshift) are the same in COSMOS as in a 
deeper observation (like HDF). Ho wever, ellipticity distribution of galaxies in a deeper survey may turn out to be 
more elliptical ()Vincent et al.ll2005l : see also the discussion on the surface brightness limit), which would increase the 
boosting factor of source ellipticity on strong lensing efficiency. We would also have a larger fraction of galaxies at the 
high redshift tail end in a galaxy survey like HDF, which increases the mean strong lensing cross section. Therefore, 
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both effects wiii produce more giant arcs. For direct comparison with observation, the observational effects, such as 
specific instrumental point-spread-function and observational seeing, are not fully included in our lensing simulations. 
As mentioned in HAA[ if the seeing was 1", then there may be only around 40% of giant arcs observed by a typical 
ground-based telescope compared with that by space-based one, such as HST. 

The same shape measurement method is used for quantifying the ellipticity of the original galaxies and lensed images. 
It could bring in some inaccuracy for measuring the ellipticity of very round images, since it is mainly designed for 
quantifying arcs. Nevertheless, this effect is rather small and negligible to our results, since we mainly focus on the 
elliptical galaxies. Besides, the intrinsic pixelization will affect the ellipticity quantification of small sources, especially 
for sources of Deff < 0.3" (see Fig. [3]). The influence of pixelization would probably shift down the measured ellipticity 
because of its relative larger effect on width measurement. Since the number of such tiny sources is small 5%) and 
the pixelization does not significantly change the ellipticity value on average (much smaller than the effect of seeing) , 
this influence would not be impor tant to our analysis. 

According to iTorri et al.l ()2004l) . cluster substructures are important for strong lensing efficiency in a major merger 
case, especially at the early state. As we adopt a softening length of 30 /i~^kpc, the softening mainly takes effect in 
the center of major merging clumps, while it would also smooth away the small substructures of a similar scale in 
a minor merger case or those remaining in a main halo. The change of the lensing cross section due to this kind of 
s moothing would not be significant in both cases. Besides, according to a co mparison plot of optical depth in Fig. 4 
of iLi et all (|2005i ). the result agrees within 25% of iWambsganss et al.l (|2004i ). in which they adopt a much smaller 
softening length of 3.2/i~^kpc than 30 h^^kpc in our case, for a source redshift of Zg — 1.0. Therefore, although the 
relatively large value of softening length could reduce the lensing efficiency, it would not affect our results significantly. 
Moreover, lensing clusters are selected from a dark matter only simulation, therefore, we cannot investigate the strong 
lensing dependence on source properties for clusters with realistic baryon distributions. 

We have only used one surface brightness limit (e.g., a detection threshold of 1.5 a) to extract the background galaxy 
images. Since the redshift of galaxy is mainly constrained by the matching with photometric catalog, it will not change 
the distribution by using a higher detection threshold in our case. A variation of detection threshold would change 
the galaxy size, but the dependence of lensing efficiency on source size is weak (see ^4.1.3p . thus the main results 
should remain unchanged. The source shape distribution also depends on the surface brightness threshold. Indeed, 
at a lower surface brightness limit, the shape of the galaxies tends to more elliptical (Vincent et a l. 2005). Thus at 
fainter surface brightness, the shape influence is expected to play a more important role on giant-arc production for 
the same lens cluster population. It is meaningful to do ray-tracing experiments by using several surface brightness 
limits to quantify the shape impact on strong lensing in different survey depths. 

The different mass range of the simulated clusters and X-ray selected ones could bring in quantitative uncertain- 
ties in the giant-arc statistics, since the strong lensing cross section strongly depend on mass, e.g., the mean cross 
section changes by nea rly 4 orders of magnitude over mass ranging from 10^^ h~^MQ to ~ 10^^ h~^MQ in Fig. 7 of 
iHennawi et al.l (|2007af l. However, the mean cross section ratio is less than 1.5 between the observed and our simulated 
clusters, according to the plot. Moreover, they use a larger mass density parameter rim.o = 0.3 and a higher matter 
power-spectrum normalization erg = 0.95 in their simulation. Besides, the simulated clusters are put at a much higher 
redshift z\ — 0.41. All these effects help to make the mass dependence of lensing probability stronger. Therefore, 
the difference of mean cross section in two comparing samples would be much less than a factor of 1.5 and becomes 
negligible in this work. Nevertheless, it is still necessary and interesting to select more comparable (between the simu- 
lated and X-ray selected) lens cluster samples in mass and redshift range for a fair comparison of giant-arc producing 
efficiency in the future. 
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